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Abstract—This paper is presenting a very efficient method for 
computing the currents on a circular loop antenna with arbitrarily 
large size. Pocklington’s integral equation is formulated for the 
circular loop geometry, and the method of moments with point 
matching is applied to cast the equation’s discrete counterpart. The 
basis functions are chosen in such a way that the relevant square 
matrix is circulant, and therefore amenable to exact eigenvalue 
analysis. Subsequently, the matrix is diagonalized and inverted an- 
alytically, yielding simple, analytical expressions for the weights of 
the basis functions, and hence the current and the input admittance 
of the antenna. The algorithm utilizes mainly elementary mathe- 
matical functions, as opposed to already known expressions in the 
literature, and yields results in the form of a rapidly convergent, 
single series, applicable to extremely large loops. 


Index Terms—Circulant matrices, current, input admittance, 
loop antennas, method of moments (MoM). 


I. INTRODUCTION 


HE circular loop antenna has been one of the most popular 
T radiators to analyze, with mathematical expressions and 
computational results presented in several pioneering papers 
[1]-[5]. It is long well-known that for very small loops the cur- 
rents is almost constant, greatly simplifying the derivation of the 
radiated field expressions [2]. As the loop radius increases, the 
current distribution becomes initially sinusoidal [6], and finally 
incorporates a large number of Fourier modes. The estimation 
of the current distribution and the input impedance/admittance 
of an arbitrarily sized loop has been the main focus of several 
works in the literature. The only way to derive any analytical 
expressions so far was the Fourier expansion of the current, 
which is clearly periodic. Invoking complex exponentials as 
basis functions in a method of moments (MoM) formulation [7] 
yields a diagonal system matrix, which is trivial to invert. Early 
papers on the topic faced convergence problems for the series 
derived. In [1], Hallen concluded that the series was divergent, 
due to a singularity occurring after the inclusion of too many 
terms. A few years later, Storer [3] avoided the large term con- 
tribution by replacing the series by an integral, which he eval- 
uated via complex analysis techniques. However, the reliability 
of the results was questioned in [4], where alternative approx- 
imations were proposed to circumvent this disagreeable situa- 
tion. Analytical expressions were derived in [4] for the weights 
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of the Fourier terms, which, however, are rather formidable, 
given only symbolically in terms of integrals of Lommel-Weber, 
Bessel and Modified Bessel functions. These expressions are 
considered classical, and have also been adopted by standard 
textbooks, such as [7]. Although exact, they are generally dif- 
ficult to compute, and are applicable in practice only for small 
to medium-sized loops. In the open literature, the largest loops 
amenable to this approach extend to six wavelengths in cir- 
cumference [8]. These elaborate expressions are still unavoid- 
able in [9], where a magnetic frill current model is used for the 
excitation. 

Alternatively, the MoM [7] can be invoked with standard 
basis and testing functions, resulting in a full square matrix, 
as opposed to the diagonal matrix required by the Fourier 
expansion. The algorithm is completely numerical, involving 
no special functions, however for very large loops the number 
of unknowns rises dramatically, rendering the method cum- 
bersome and inefficient. Fast methods [10] can obviously 
be employed to accelerate the MoM algorithm, however the 
unknown current is still computed through a numerical matrix 
inversion, since no closed form expressions can be derived. 
Therefore, ill-conditioning may deplete the methods’ efficiency 
for very large loops, since the number of iterations required for 
the linear system solution may increase significantly. Conver- 
gence checks are thus laborious, whereas loop design becomes 
a very tedious task, since several time-consuming repetitions of 
the algorithm may be required to obtain the desirable antenna 
parameters. 

All these complications could be resolved if simple, analytic 
or semi-analytic expressions could somehow be derived for the 
current and the input impedance of the antenna. In this paper 
it is shown that such expressions can be extracted indeed, by 
invoking the properties of circulant matrices [11]. This type of 
matrices can be diagonalized and inverted in closed form, irre- 
spective of their dimension. The important implications of this 
property can be readily seen in scattering analysis of structures 
with circular symmetry [12]-[14]. In particular, [12] elaborates 
on the use of block-circulant matrices for the radiation analysis 
of circular arrays of dipoles, [13] utilizes similar concepts for 
the scattering analysis of discrete bodies of revolution, whereas 
[14] is a representative of a series of papers by the same research 
group on the optimization of the Method of Auxiliary Sources 
applied to circular cylinders of various materials. The geometry 
of the circular loop is clearly amenable to such an approach, pro- 
vided that the basis and testing functions of the MoM formula- 
tion are suitably chosen. It will be proven that both the current 
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Fig. 1. 


Geometry of a thin circular loop antenna. 


and the input admittance of the loop can be given in terms of 
a rapidly convergent, single series, which involves almost ex- 
clusively elementary functions. The advantages of such expres- 
sions cannot be overemphasized, since all aforementioned draw- 
backs of other methods are completely circumvented. Accurate 
and robust calculation of the antenna parameters for loop radii 
of hundreds of wavelengths becomes feasible in a few seconds, 
even on a scientific pocket calculator. Furthermore, convergence 
checks are trivial to perform, leading to high precision of the nu- 
merical values obtained. 

The structure of the paper is as follows: Section II formu- 
lates the MoM system in a suitable way, the matrix is inverted 
exactly, and semi-analytical expressions for the current and the 
input admittance are derived. Section III presents several numer- 
ical results for loop circumference up to 100 wavelenths, while 
comparisons with reference data for small loops validate the al- 
gorithm and the computational code. Section IV outlines and 
concludes the paper. 

A eJ¥t time dependence is assumed and suppressed 
throughout the paper. 


II. ANALYTICAL FORMULATION 
A. Matrix Fill 


Assume a thin circular loop of radius 6 and wire radius a << 
b, driven, without loss of generality, by a delta gap generator 
(Fig. 1). Let the unknown current on the loop be I(¢), which is 
assumed to flow on the wire axis of radius b, given that the wire 
is very thin. The azimuthal component FE’ (r) of the radiated 
electric field at an arbitrary observation point r(r,@) is given 
by [7] 


Qn 
r ee / / , Coat! i 
Ex(r) = —jwpr iC ) cos(¢ — ears: 
0 
20 
| 1 dI(¢’) d e dk lr—r" | i 
" iwer dd! dd ae dp (1) 
0 


where, as usual, |r — r’| is the distance between source and ob- 
servation points, w is the circular frequency, ¢ is the dielectric 
permittivity of the medium, jz is the magnetic permeability and 
k; is the wavenumber. Applying the appropriate boundary con- 
dition for the electric field at a discrete set of N points (point 
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Fig. 2. The step-pulse basis function, which adequately approximates a 
triangular basis function (dashed line). 


matching) with azimuthal positions $m = (2r/N)m = m@, 


m =1,...,N on the outer surface of the wire yields 
20 : 

i . e—skR 
(6-44) ES (bm) = jhniXb-+a) | T(§) c08(bm 8) Fa 
0 

af do") e~H#® 
e 
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where Ei(¢m) is the azimuthal component of the incident elec- 
tric field due to the antenna excitation calculated at ¢,, 


R= bs (*5*) +4 (7) (3) 


and 7 is the medium intrinsic impedance. 

To complete the MoM formulation, and cast the linear system 
of equations, the current J(#) must be expanded into a weighted 
sum of basis functions, i.e. 


N 
1(¢) = >> InFu’). (4) 
n=1 


It should be kept in mind that the purpose in our procedure is to 
obtain a circulant matrix, therefore the basis functions F;,(¢’) 
should be chosen accordingly. Pulse basis functions would be 
suitable from this aspect, but the accuracy of the final result 
would be rather poor [7]. Triangular basis functions would guar- 
antee higher accuracy, however they would unnecessarily com- 
plicate the matrix elements calculations. A very good compro- 
mise between accuracy and simplicity is the step-pulse basis 
functions [7], which adequately approximate a triangular basis 
function, yet maintaining the salient features of a regular pulse 
(see Fig. 2). The analytical expression for the step pulse of Fig. 2 
is 


Fu(#!) = Tn($') + ST n-4O) + ST) ) 
where 


m,(6!) = 41 
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From (2) and (4) it becomes evident that the derivatives of 
F,,(¢’) must also be utilized in the formulation of the linear 
system. Using the properties of generalized functions (distribu- 
tions), (5) and (6) yield 


gp -to(v (0-93) “9(0-(0-4)3) 
“H(e-(o+4)3)-Be #38) 0 


where 6(e) is the Dirac’s delta function. Again, using the prop- 
erties of the delta function, and approximating the 0/0¢ deriva- 
tives by finite differences, the second integral in (2) is approx- 
imately equal, due to (4), to a weighted sum of the following 
expressions: 


Qn . 
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Ray = b asin? VE 5 (2)? (9) 
KY — ») b 


From (8) and (9) it is clear that this portion of the MoM ma- 
trix elements, corresponding to the second integral in (2) is a 
N-periodic function of the difference m — n only, and there- 
fore demonstrates circulant behavior, as required. Similarly, the 
first integral in (2) can be approximately calculated. In the case 
where m # n, the evaluation is straightforward, yielding 
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J Fa(#))008(m - 
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which is again a function of m — n only. Form = na much 
more refined approach is required, due to the behavior of the 
Green’s function kernel, and the resulting very small values of 
the denominator. The two last terms in the right-hand side of 
(10), which correspond to the side pulses in Fig. 2 are still valid, 
however the first term, which corresponds to the central pulse, 
must be evaluated with higher precision. Local linearization of 
the wire, as proposed in [7], turns out to yield inaccurate results. 
On the contrary, a more careful integration in polar coordinates 
yields ({15, p. 173, 2.597,1] 
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where F is expressed as in (3), evaluated at 6 = @m, and 
_ VW1472sin ¢ 2b 
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pd 
a 
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(elliptic integral of first kind). (12) 


The elliptic integral encountered in (12) is the only nonelemen- 
tary function involved in the entire formulation, however it is 
very easily computable through efficient and simple subrou- 
tines, readily available in the literature, e.g. [16, p. 260]. 

In brief, the MoM linear system produced by (2) is written on 
the basis of (8)—(12) as 


{V} = [4] {2} 
where {I} is the vector of unknown weights, V, = 
(b+ aE (nd) and [Z] is a square matrix, whose elements 
are N neriedic functions of the index difference p = m—n 
only. This is the defining property of a circulant matrix, namely 


(13) 
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B. Analytical Solution of the Linear System 


Since [Z] is circulant, its eigenvalues A, and eigenvectors 
{g}q are given by the theorem in Appendix A, therefore [Z] can 
be explicitly diagonalized, and therefore analytically inverted. 
Elementary algebraic manipulation, similar to [14], based on the 
expressions of Appendix A, yields the exact solution of (13) 


ened Qa 
In = HDL Fo {ulema D | Vo 
m=1,...,N (17) 


for arbitrary excitation {V}. For a delta source of voltage V 
across the gap, located at ¢ = 0, evidently 


6n,1 


(18) 
where 6,,1 is the Kronecker delta, hence (17) is simplified as 


‘a a ~ 1 20 
Im = 95 (1 3 5) 2a, exp { tm = 1N(q- D5} (19) 


and the input admittance is therefore given by 


N 
Yn pag (4 UG, 
2 


The physical meaning of (20) is that the loop antenna can be 
characterized as a circuit of parallel impedances, equal to the 
eigenvalues of the MoM system (13) within a constant. Expres- 
sions (17), (19), and (20) constitute a great simplification with 
respect to standard formulas for the current and the input ad- 
mittance given in [3], [4], and [7], since almost all calculations 
can be performed through single sums of elementary functions. 
It is pointed out again, as in the introduction, that the respective 
expressions in [3], [4], and [7] utilize complicated integrals of 
special functions that are cumbersome to compute, especially 
for large loops. Admittedly, the method proposed herein is not 
purely analytical, in a sense that the resulting expressions are 


es, 1 V a\ VN 
Vn = (b+ a) E4(nb)bn,1 = (b+4) 55 4n,1 = 1+—-)—— 
2 ( 5) Qn 
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Current amplitude convergence, kb=1 


current amplitude (A) 
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Fig. 3. Convergence of the current amplitude for a loop with kb = 1 and 
b/a = 15. 


neither “exact,” nor “elegant,” such as the ones in [3], [4], and 
[7]. However it will be shown in the next section that, from a 
computational point of view, (19) and (20) are much more effi- 
cient and highly accurate. 

Furthermore, the capability of analytical matrix inversion, 
and the consequent insensitivity to matrix conditioning, offers 
to the method several unique advantages over purely numer- 
ical algorithms. Indeed, unlike the formulation proposed in this 
paper, ill-conditioning may decrease the accuracy of numer- 
ical methods combined with direct solvers, and reduce the con- 
vergence rates of methods based on iterative solvers. In terms 
of memory requirements, the demands are obviously very low 
(O(N)). Finally, in terms of computational cost, the entire pro- 
cedure is at least comparable to fast methods, even in the ideal 
case for the latter, i.e. when a single iteration is sufficient. This 
is evident from (A1), (17), and (19), which essentially represent 
discrete Fourier transforms, and which can rapidly be calculated 
through FFTs, thus reducing the complexity of the algorithm 
from brute-force O(N”) down to O( Nlog2N). Outperforming 
fast methods, however, the result is given in only one step. 


II. NUMERICAL RESULTS 


To validate the analysis presented in Section II, several nu- 
merical results are extracted and compared to reference data, 
whenever available. First of all, the convergence rate of the se- 
ries expressions involved is examined. Fig. 3 depicts the conver- 
gence behavior for the current amplitude on a loop with kb = 1 
(ie. with circumference equal to A) and b/a = 15. The con- 
vergence behavior for the current phase is shown in Fig. 4. In 
both figures the horizontal axis corresponds to the azimuth angle 
along the loop circumference. From both figures it is concluded 
that convergence is achieved for a number of unknowns approx- 
imately equal to N = 40, i.e. for an element arc length equal to 
\/40. Of course, the radiation patterns are not as sensitive to the 
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Current phase convergence, kb=1 
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Fig. 4. Convergence of the current phase for a loop with kb = 1 and b/a = 
15. 
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Fig. 5. Current amplitude for loops with various circumferences and b/a = 
23.62. Compare with [18, p. 185, Fig. 5.7]. 


element arc length as the actual currents [17], and they usually 
converge must faster, typically for element arc lengths equal to 
\/10. 

Comparison with reference results is performed by applying 
the analysis to loops of increasing circumference (ranging from 
0.1 to 0.4 wavelengths) and radii ratio equal to b/a = 23.62 
(or 21In(27b/a) = 10, to comply with the notation used in [3] 
and [18]). The results for the amplitude are shown in Fig. 5 and 
for the phase in Fig. 6, and are clearly identical to [18, p. 185, 
Fig. 5.7], borrowed in turn from [3]. Moreover, results for the 
input admittance are depicted in Figs. 7 (conductance) and 8 
(susceptance). Fig. 7 is identical to [7, p. 93, Fig. 5.5], which 
was plotted using the standard, cumbersome expressions. Fig. 8 
is very similar to [7, p. 93, Fig. 5.6], although not identical. The 
rather insignificant discrepancy is fully explained by the fact that 
the input susceptance for a delta gap model cannot be rigorously 
defined, since it approaches infinity as the gap width vanishes 


Current phase (rac) 


Oo 20 40 60 80 100 120 140 160 180 
azimuth angle phi (deg) 


Fig.6. Current phase for loops with various circumferences and b/a = 23.62. 
Compare with [18, p. 185, Fig. 5.7]. 
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loop radius b {in wavelengths) 
Fig. 7. Input conductance as a function of the loop radius, for various ratios 


b/a. Compare with [7, p. 93, Fig. 5.5]. 


[5]. Therefore, any susceptance calculations strictly do not con- 
verge [7], and are meaningful only as approximate estimates for 
realistic loops [5] with excitations modeled by small, but non- 
vanishing, gaps. 

Since all comparisons with reference data showed excellent 
agreement, validating the analysis, the semi-analytical results 
in Section II will now be used for current computations of very 
large loops. Figs. 9 and 10, depict the current amplitude and 
phase respectively, as a function of the azimuth angle along the 
loop circumference, for a = 4/100 and kb =10, 25, 50, and 
100. Convergence was carefully examined, being achieved for 
N = 600, 2500, 5000, and 8000, respectively. Finally, Fig. 11 
depicts the dependence of the input conductance as a function of 
the loop circumference, the latter extending to very large loops. 
Input susceptance is not plotted for reasons explained earlier. 
All calculations, including convergence checks, were performed 
on a typical PC in a matter of a few seconds, based on an ele- 
mentary computer code. Hence, the method is valuable for very 
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b/a. Compare with [7, p. 93, Fig. 5.6]. 
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Fig. 9. Current amplitude as a function of the azimuth angle for loops with 
various, very large circumferences and a = \/100. 


large loop designs, since it is very simple to formulate, and since 
it retains the accuracy of the purely numerical MoM, requiring, 
though, very low memory and CPU time resources, and also 
avoiding instability problems caused by possible matrix ill-con- 
ditioning. The results shown in Figs. 9-11 are related to loop 
sizes far larger than the 6. previous upper bound, recently doc- 
umented in the literature [8]. 


IV. SUMMARY AND CONCLUSIONS 


A semi-analytical formulation for the efficient and accurate 
calculation of the currents on extremely large loops was pre- 
sented. The basis of the analysis was the MoM, whose linear 
system was appropriately cast in such a way that the relevant 
matrix was of circulant form. The matrix was therefore diag- 
onalized and inverted by an exact procedure, yielding analyt- 
ical expressions for the currents and the input admittance. The 
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Fig. 10. Current phase as a function of the azimuth angle for loops with 
various, very large circumferences anda = \/100. 
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Fig. 11. Input conductance as a function of the loop circumference, extending 


to very large loops and a = \/100. 


formulation utilized only elementary mathematical concepts, as 
opposed to standard analytical methods, previously developed 
in the literature, which involve complicated integrals of spe- 
cial functions. Moreover, the analytical invertibility of the MoM 
matrix reduces tremendously the MoM computational cost, with 
respect to memory and CPU time requirements, a very useful 
feature for the analysis and design of very large loops, where 
repetitive convergence checks are computationally very expen- 
sive and tedious. Furthermore, as opposed to fast methods, the 
final results cannot possibly be affected by ill-conditioning of 
the MoM matrix and hence, are always numerically stable. The 
method is not restricted by any upper bound of the size of the 
loop, and accurate data for circumferences up to 100 wave- 
lengths were plotted, i.e. for loop sizes much larger than the 6 
wavelength limit documented in the literature so far. 
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APPENDIX A 
Theorem [11]: Let[A] bea N x N circulant complex matrix, 
ie. [A] = circ(a,,a2,...,an). Then [A] is normal, and thus 
diagonalizable. Its eigenvalues are given by 
N fs 
vq = So apwP-DE-Y, q=1,...,N,w=en 
p=1 


(Al) 


and the corresponding eigenvectors by 


[1,00-0, wD, jut —DO-D) (a2) 


= ok 
{8ha = Ty 
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